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Abstract. The average density of zeros for monic generalized polynomials, P n {z) = 4>{z) + Ylk-i c kfk(z), with 
real holomorphic <f>, fk and real Gaussian coefficients is expressed in terms of correlation functions of the values of 
the polynomial and its derivative. We obtain compact expressions for both the regular component (generated by the 
complex roots) and the singular one (real roots) of the average density of roots. The density of the regular component 
goes to zero in the vicinity of the real axis like |Imz|. We present the low and high disorder asymptotic behaviors. 
Then we particularize to the large n limit of the average density of complex roots of monic algebraic polynomials of 
the form P n (z) — z n + 2fc=i c kZ n ~ k with real independent, identically distributed Gaussian coefficients having zero 
mean and dispersion S = The average density tends to a simple, universal function of £ = 2nlog \ z\ and A in the 

domain £ coth | <C n\ sinarg(z)| where nearly all the roots are located for large n. 



> 
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I. INTRODUCTION 



o 

\Q • Let P n be a (monic) algebraic polynomial of degree n, 

o , , 

\D . Pn(z) = z n + a 1 z n - 1 + --- + a n . (1.1) 

The roots of P n are (algebraic) functions of its coefficients a\, d2, • • • , a n . If the coefficients are random variables, 
then the roots of the polynomial will also be random variables. The study of the distribution of the roots of random 
polynomials began with the investigation by Bloch and Polya continued by Littlewood and Offord Kac Q] 
and many others (for reviews see |)|]6|]) of the number and distribution law of the real roots of random polynomials 
with real coefficients. 

Most investigations of the distribution of zeros for random algebraic polynomials have dealt either with the initial 
problem - the real zeros of real polynomials - or with the complex zeros of complex polynomials. In the latter case, 
one is helped by the fact that it is possible to transform the joint distribution function of the coefficients into the joint 
• i-H , distribution function of the roots explicitly using the relations between the coefficients of an algebraic polynomial and 
' its roots.[] The average distribution of roots for homogeneous algebraic polynomials, 

| n-l 

Pn-l(z) = ^CjZ J , 
3=0 

with real normal Gaussian coefficients Cj, was recently investigated by Shepp and Vanderbei^j ]To| . They obtain a 
generalized Kac- Rice formula for the density of complex roots and explore its large n limit. The fraction of the 
expected number of roots contained in an angular sector S{0\,02) which does not intersect the real axis^j tends to 
I #2 — 0x|/2?r. Most roots are concentrated in a small annulus (with width ~ n _1 ) near the unit circle. Af(R), the 
expected number of roots in a ball of radius R, satisfies 



1 Which follow by comparing the polynomial P n (z) = ~Y^l =(l o,kZ n ~ k with its roots expansion P n {z) = aon^=i( 2 ~ z i)- The 
Jacobi determinant for this transformation was computed by Girschick and Hammersley 
2 We thank an anonymous referee for bringing this paper to our attention 
3 The analogous result for complex coefficients was obtained by Shparo and Shur O], 
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lim -W(e^) = — - -. (1.2) 

n— >oo n 1 — e s s 

In a recent preprint by Ibragimov and Zeitouni this result is generalized to coefficient distributions belonging to the 
domain of attraction of a a-stable law ]l^] . 

In this paper we study the distribution of the roots for more general random polynomials with real coefficients. The 
initial motivation for the present work was given by some interesting properties of the roots of the Szego orthogonal 
polynomials jl3| related to the Wiener transfer function [E3. Szego polynomials are orthogonal polynomials in 
the non negative measure on the unit circle of the complex plane d/i(#) generated by a non-negative Toplitz form 
T{k - I) ;k,l = 0,1, ■■■■ Here 



T(k) = r d^ey 

J —7C 



e- Ifctf (1.3) 

are the measure's moments. In the case of the Wiener transfer function the moments are equal to the autocorrelation 
function of a discrete finite real signal sample (X(p) — 0; p < or p > N): 

N-l 

T{k)=F N {k)=Y J X{p)X{p+k) (1.4) 

Let the signal be the sum of a useful signal, consisting of several harmonic components with frequencies oj v , and noise. 
Then there is strong numerical and some analytic evidence that for signal-to-noise ratios that are not too low, the 
roots of the Szego polynomials of order n with 1 <C n >C N break into two groups. The first one converges rapidly to 
the unit circle at the points e lw ", [p^LpTf. The other group, which also converges to the unit circle, but more slowly, 
is nearly equispaced - resembling an one-dimensional crystal. It has universal statistical properties even when the 
harmonic components of the signal are absent ||||. 



If the useful signal is absent, the Szego polynomial of order n is well approximated by (1.1) with at — N 2 c k with 



approximately Gaussian coefficients Ck if 1 <C n -C N [18|. In the presence of a useful signal, the polynomials have 
a mo re complicated form, which can be obtained by substituting more general polynomials for the monomials z k in 
(1.1). Thus, we are led to considering generalized random polynomials. 



Let 

<f>{z), fk(z), fe = l,2,...,n , neN (1.5) 

be holomorphic and linearly independent functions of the complex variable z in a domain of the complex plane, 
z = x + iyEDcC. Let c(w) be a random n-dimensional vector with components c k {uj), k = 1, 2, . . . , n , n G N. 
Here lu G Q where Q is a probability space. 

We define random generalized monic polynomials of degree n by 

P„( 2 ) = 0(z)+ 2 c k f k (z). (1.6) 

l<fc<n 

Homogeneous polynomials (of degree n — 1) correspond to setting cf> = 0. Alternatively, a monic polynomial of 
degree n can be regarded as a homogeneous one, having a singular distribution for the coefficient c„+i {8{c n +i — 1) ). 



Setting <j){z) = z n ,fk{z) = z n , fc = 1, 2, . . . , n , n G N we obtain the random monic algebraic polynomials, (1.1) 



Trigonometric, hyperbolic and other types of random polynomials p[ can be obtained by su itab ly defining and f k 



In the following we will often omit the qualifier generalized and call the objects defined by Eq.(L6) simply polynomials 



Let us note that any non-random affine transformation of the random vector of the coefficients, 

Cfc — > Qfc + Kfc TO c TO (1.7) 
wher e a is a constant vector and K a non singular matrix, transforms the polynomial P n into one of the same form 



(1.6) but with the basis set (1.5) replaced by 



m rn 

Thus, by redefining the basis set, we can always consider that the mathematical expectation of the polynomial is 
equal to its deterministic part, <fi. 
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Let F(z;u>) be a random holomorphic function of z € D C C, i.e. a family of functions F(z;u>) indexed by u> G 51 
which are almost surely holomorphic for z in the domain D. Let 



z r {bj) = x r {oj) + iy r (u>), r 



(1.9) 



be the solutions of the equation 



F(z;w) = 0. 



(1.10) 



The zeros of F are random variables. In each compact subdomain D± C D, there may be only a finite number 
N{D±;uj), of zeros for each realization u> G f2, since the accumulation points of zeros cannot lie inside the domain of 
holomorphy. 

The density of zeros of the random function F is the random distribution on M 2 : 



p(x, y;uj) = ^2 s ( x ~ x r(u))8(y - y r {^)) 

r 



(l.n) 



where S and 8 ^ are, respectively, the Dirac distribution on R and on K 2 . Using the definition of the Dirac 6 we may 
rewrite ( |l.ll| ) as 



p(x,y;u) = F(z;lu) 6™ (F(z; to)) , 



where the Cauchy-Riemann conditions were used to calculate the Jacobian. The compact notation 

for the derivative of F will be used throughout this paper. 

The expected (average) density of zeros of F is obtained by averaging over the realizations: 

T>(x,y) = E{p(x,y;o)} . 

In a similar way, one may define the two-point correlation function of the zeros by 

^2(^1,2/1; ^2, 2/2) = E{p(x 1 ,y 1 ;°)p(x 2 ,y2;o)} 



(1.12) 



(1.13) 



(1.14) 



(1.15) 



and higher, m-p oint, correl ation functions for the roots. 

Substituting (1.12) into (1.14) and introducing the joint distribution function of the values of the function and its 
derivative at the point z = x + ly 

V(a,a;z) = E |(5 (2) (a - F(z; ))S (2) (a - F{z; o)) } , (1.16) 

we see that the average density of roots at the point z = x + \y is given by the Kac-Rice [p|j5|pct| formula 

V(x,y)=E{\F(z;o)\ 2 8( 2 \F(z;o))} 

|a| 2 7>(0,a;z)d (2) 5. (1.17) 

in terms of V(a,a; z). Similar formulas, involving only the joint distribution function of the values of the function 
and i ts derivative at the selected points may be written fo r the m-point correlation functions of the roots such as 
( 1.15 ). We will s tudy the average distribution of roots, Eq.( 1.14 ), in the case when F(z,uj) is a generalized random 
polynomial, (1.6), with Gaussia n co efficients. 

If the components of c(uj) in ( JLg ) are Gaussian (real or compl ex) ra ndom variables, the joint distribution function 
of the values of the polynomial and its derivative at some point z (1.16) will also be also a Gaussian distribution which 
is determined by the correlation functions (covariance matrix) of the values of the polynomial and its derivative at 
that point. This means that we may compute explicitly the integral in Eq.(1.17) for the average density of roots in 
the Gaussian case and the corresponding expressions for the m-point correlation functions of the roots. 
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In Section [II we will obtai n th e average density of roots, Eq.(1.14), for real Gaussian generalized monic polynomials, 
when the basis functions in (1.5) are of real type 



[4>{z)Y = d>(z% [h{z)T = f k (z*). 



(1.18) 



Here * denotes complex conjugation. The coefficients Ck are independent, identically distributed ((iid)) real Gaussian 
random variables with zero expectation value. The assumption of a joint Gaussian distribu tion is important, while 
the restriction to the iid case is inconsequential. Indeed, by a suitable linear mapping of type (1.7) any finite Gaussian 
distribution may be mappe d on to the standard one {iid with zero average and unit dispersion). By the above remark 
this leads to a redefinition (l.S) of the basis set. 

We will obtain a ge neral formula giving the expected density of complex roots at points with Im(z) ^ 0. Due to 
the reality condition, ( [1.18] ), there is also a singular component of the expected density of roots located on the real 
axis. We will obtain a generalization of the Kac-Rice formula for it. For homogeneous polynomials with non-central 
Gaussian distributed coefficients we recover a result by Edelman and Kostlan [6], who also pointed out that taking 
a singular correlation matrix limit for that case will yield the density of real roots in the monic case. We will show 
that the density of complex roots tends to zero, like |Im(.z)|, in the vicinity of the real axis. In the high randomness 
limit (large dispersion S of the Gaussian distribution) the expected density of roots for monic polynomials approaches, 
obviously, that for the homogeneous ones (4> = 0). In the weak randomness limit, when the dispersion S goes to zero, 
the roots concentrate near the zeros of the deterministic part cf>. Near simple zeros of 4> the distribution tends to a 
Gaussian one. Near zeros of <p with multiplicity k > 1, it has generically 2k maxima at a distance ~ 5t° from the 
zero's position. An interesting situation arises when the multiplicity of the zero is large. 

Since only the joint distribution function of the polynomial and its derivative is assumed to be Gaussian, the results 
of Section III could also be relevant to cases when the coefficients are non-Gaussian but one may prove a limit theorem 
for the joint distribution function. 



In Section IV we will apply the general results of Section III to the density of roots of monic algebraic random 
polynomials, (IT), with iid coefficients having a real Gaussian distribution with zero average and dispersion 5 = 
in the large n and N limit. In this case the deterministic part has a real zero of multiplicity n at the origin. As 
mentioned above, in |l8[ the Szego polynomials associated to the Wiener transfer function for pure noise were shown 
to have this form if the sample length is N n ^>> 1 . 

Defining a new rescaled coordinate £ = 2nln|z| and rescaling also the dispersion, S^ 1 = N = n\, we obtain a 
simple asymptotic expression for the average density of roots valid for large n at points which satisfy 



£ coth - <C n\ sinarg(z) 



(1.19) 



This is a small neighborhood of the unit circle which does not intersect the real axis and it contains nearly all the 
zeros. To leading order, the average density in this domain does not depend on arg(z). It approaches n 2 times a 
universal function of £ and A. For sin# ^= 0, 



lim -^T>(e^ 




1 



exp 



-A- 



(1.20) 



In the case of homogeneous polynomials, which corresponds to A = 0, ( 1.2C ) is equivalent to (F2), the result of 
Shcpp and Vanderbei |lC| ]. 

We will finally remark that the asymptotic distribution of complex roots for polynomials with real Gaussian coeffi- 
cients, (1.20), which do not have the rotational symmetry of the polynomials with complex coefficients, nevertheless 
coincides with the rotationally invariant one for the latter, which can be readily obtained using the methods we use. 
For the complex homogeneous case (A = 0), the asymptotic estimate (F2) is due to Arnold 21 |]. 

A part of the calculations for the algebraic polynomials will be presented in more detail in the Appendix. Before 
proceeding further, we will defin e in the next section some notations that will allow us to write in a compact way the 
joint distribution function (1.16) and our results when the random function is a polynomial with Gaussian coefficients. 



Acknowledgments GAM, DB and JDF acknowledge the support of the European Science Foundation for par- 
ticipation at the Como Workshop Classical Mechanics Methods in Quantum Mechanics; GAM and GM thank Alfred 
Msezane and Carlos Handy for hospitality at CTSPS, Clark Atlanta University, where parts of this work were done. 
JDF and GM acknowledge support from the Programme Galilee de collaboration scientifique France - Italic 
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II. A DIGRESSION ON NOTATIONS 



For any holomorphic function, h, we will write 



for its derivative. The real and imaginary parts of a complex number or function h will be denoted by the subscripts 
1 and 2 respectively. 

h x = Re(h), h 2 = lm(h). (2.2) 
As remarked in the Introduction, for a given realization the coefficients define a point in R™, while the basis set of 



functions, </> and k = 1 . . . n, maps C ~ R 2 in C x C" ~ R 2 x R 2 ™. The values of the polynomial P n (z), Eq.(1.6) 
and its derivative P n {z) are affine mappings (with complex coefficients) of the real n- vector of the coefficients into 
C ~ R 2 . Let us introduce more compact notations for the various vector structures we are dealing with. 

We will use the same lower case Greek letter (e.g. <p) for both the complex scalar 4> S C and the real 2-vector <j> 
whose components (<f>x, <f>2 ) £ R 2 are the real and imaginary parts of the complex scalar <f>. 

Boldface letters denote real n-vectors R n 3 f = (fx, fi, ■ ■ ■ , f n )- When it does not lead to ambiguities, we might 
abuse the notation to denote a complex n- vector - like the vector made of the basis set of functions fk(z), k = 
l,2...,n - by f(z). In such cases, the real and imaginary parts are f\(x,y), f2(x,y) and we may write f(z) = 
fi(x,y) +ii 2 (x,y). 

Latin letters, /, g, will be used for real 2rt-vectors with components /fe Q , gk a , k = 1, 2, . . . , n, a = 1,2. We use the 
same letter in the typefaces mentioned above for related objects, like f(z), (fx(x,y), f 2 (x,y)), f(x,y) for the basis 
set of functions at point z considered as a complex n-vector made of two real n-vectors or as a real 2n- vector. 

Greek subscripts always run over the set {1,2} while the Latin ones over {1, 2, . . . , n}. 

For a Ga ussi an distribution of the coefficients, the joint distribution function of the values of the random polynomial 



P n (z), Eq.(1.6), and its derivative P n (z) is determined by their correlation matrix. The calculation of averages can be 
done using Wick's theorem, replacing the products of coefficients by their expectation value. This reduces further to 
a contraction over the Latin indices if the Gaussian process is a direct sum of normal ones.^ The compact notations 
defined in the rest of this section will help us with the necessary bookkeeping. 

We use bold square brackets [o, o] for the contractions over the Greek indices and bold round brackets (o, o) for the 
contractions over the Latin ones. Using these conventions, 

2 

[&V] = X>°^*' (2-3) 

a=l 

will be the real scalar product of the 2-vectors <fi and ip] [</>, /] will denote the real n-vector with components 

2 

a=l 

obtained by contacting the direct product between the 2-vector <j> and the 2n- vector / over the Greek indices. In 
an analogous way, [/, g\ will be the second order tensor obtained by contracting over the Greek indices the direct 
product of the 2n-vectors / and g. Its components form the n x n matrix 

2 

[f>9] jk = 1 ^2fja>gko t , (2.5) 

a=l 

In a similar way, the scalar product between the real n-vectors f and g is 

n 

(f,g)=X>0k, (2-6) 

k=i 



As noted in the introduction this is always possible choosing a suitable (1.7) 
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while 

n 

(f,g) al3 = ^2fk a gk/3, (2.7) 
fc=i 

are the elements of the 2x2 (real) matrix (/, g) , while 

n 



k=l 



are the components of a 2-vector. 

Finally, we use matrix notations like (f,g)ip for the product of the 2x2 matrix (/, g) and the 2-vector ip, 

2 

((/^) Q = E(/'SU#- (2-9) 

,3=0 



III. AVERAGE DENSITY OF ROOTS FOR RANDOM GAUSSIAN GENERALIZED MONIC 

POLYNOMIALS 

In this section we consider monic holomorphic polynomials of type Jl.of ) with coefficients distributed according to 
a Gaussian law. As we remarked in the Introduction, the case of homogeneous (non monic) polynomials may be 
obtained by setting (j> = 0. There we also noted that it suffices to consider only the case of iid normal Gaussian 
coefficients (with zero expectation and unit dispersion): 

E{c fc } = 0, 

j, A: = 1,2, ...,n. (3.1) 

E{cfeCj} = 6 k j, 

The calculation of the average density of roots, their two-point correlation function T>, T>2 and even higher correlation 
functions T> m may be done in closed form since all the integrals will be Gaussian, although the formulae will tend to 
become rather encumbering with increasing m. 

A. Average density of complex roots 

The starting point will be the formula ( 1.17| ) for the average density of roots for the polynomial 



P(z) = <f>(z)+F{z), (3.2) 
where (j) is the deterministic (non random) part and 

n 

F(*)=£ Cfe / fe (2) = (c, /(*)). 
fe=i 



Here c is the real vector of the Gaussian coefficients, and we use the bold round bracket notation (2.8) introduced in 
the previous section. 

For the sake of clarity of the main points of the calculation, we will first illustrate our approach in the case of 
homogeneous (non monic) polynomials, i.e. 

= 0, 



obtaining a generalization of the results derived by Shepp and Vanderbei |10| in the algebraic case {fj{z) = z J_1 ). 
Subsequently we will obtain the average density of complex roots in general case of monic polynomials. 

In the homogeneous case the joint distribution function for the polynomial and its derivative, P and P, at the point 
z = x + iy coincides with the one for F and F at the same point. Then, 



G 



(3.3) 



Here we used the bold square bracket notation ( |2.3| ) for the contractions introduced in section || to rewrite \a\ 2 as 
[a, a]. 

Now, V will be a Gaussian distribution determined by A(z), the 4x4 correlation matrix of F(z) and F(z). 

( • • ' " 

Its elements are the expectation values of products of elements of the 4-vector I Re(F), Im(F), Re(F), Im(F) 

(Fi ,F 2 ,Fi,F^j . For example: 

A M (z) = E{Re(F(z)) Im(F(z))} = £ Ai(z) /,- 2 (z)E{c feCi }. 



Using (3.1) this reduces to 



A 14 (z) = X)/wW) Aa(z). 



We may use the bold round bracket notation (fO|) to write the 4x4 matrix A(z) as a 2x2 block matrix whose elements 
are 2x2 matrices 



A(z) 



U(z),m) (m,m) 

(f(z)J(z)) (f(z)J(z)) 



(3.4) 



As a correlation matrix A(z) is symmetric and non-negative. Th e sam e is true for the diagonal blocks in fl3.4| ). 
Since the functions fk{z) are linearly independent and of real type, (1.18), A(z) will be generically non singular off 
the real axis. For Im(z) ^ 0, we may then write the joint distribution function of F and F 



V{a, a; z) 



(27r)V det A (^) 



exp 



(3.5) 



a 2 



Q'2 



Here a = [ ai I and a — ( < ? 1 I are 2-component vectors and we use T to denote the transposed. 



For Im(z) ^ 0, where the rank of A is equal to four, we may introduce the block-Cholesky decomposition of the 
positive matrix A into block-triangular 2x2 factors: 



A = K K, 



where the upper block triangular matrix 



K = 



(/,/)' (/,/) 5 (/,/)' 

(A/). 



Here we introduced the notation 

(A f) ± = (/ - T ± f, f - T ± f) = (/, /) - (/, /) (/, fY\f, f) , 

with the matrix 

The square root of a matrix, M , is defined by the well known relation 



(3.6) 

(3.7) 

(3.8) 
(3.9) 



M5 = - 

7T 



dt 



l + i 2 M-^ 



which is true as long as M has no eigenvalues on the interval (— oo,0]. For brevity's sake we also dropped the 
dependence on z = x + ly, which we will continue to do whenever this does not lead to ambiguities. 
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If the matrix A is strictly positive, its diagonal block (/, /) has the same property. Then, 



where 



K = 



A" 1 = K-^K- 1 )' 



(/,/) 5 -(/,/) *(/,/) (f,f) ± 



(/>/). 



(3.10) 



The matrix (/, /) was defined above, (ft.S| ). 



Substituting ( |3.5|) into (3.3) we may rewrite (^3) as 

1 



V 



[a, a] 



exp 



(2n) 2 VdetA 

This Gaussian integral is readily evaluated yielding 



((/,/) J 



/det(/,/) ± 



(X 



(3.11) 



But from (3J3) and (3/7) and taking into account the block-Cholesky decomposition introduced above 

det A = det K 2 = det (/, /) det (/, f) ± . 

Thus, we get a simple formula for the average density of complex roots extending the Shepp and Vanderbei one |Tc| l 
to random Gaussian homogeneous generalized polynomials: 



V 



Tr(/,/) 



2tt Jdet (/, /) 



(3.12) 



Let us now return to the general case of monic polynomials, when the deterministic part <j) is not identically zero. 
The joint distribution function of the values of the polynomial and its derivative P and P coincides with the one for 
F and F shifted by their deterministic parts <j> and, respectively, <fi. Thus, for Im(z) ^ 0, 



V = J [a,&]V(-(f>,&- <j))d {2) &, 



where the function V is again given by ( |3.5| ): 

1 



(27r) 2 v / detA 



exp 



4< 



6 T ,a T -0 T )A" 1 



a — <p 



(3.13) 



(3.14) 



The inhomogeneous Gaussian integral ( 3.13] ) with V given by (3.14) is calculated in a similar way to the previous one 
by using the block-Cholesky decomposition (3.7-3.10) and introducing a new integration variable 



= a - <f>±, 

where 

<f>± = 4> - T_|_<^. 
Finally, the expected density of roots in the monic case is 



V 



exp {-± [<£,(/,/) V]} 



2^^/det (/,/) 

where (/,/), and Tj^ were defined above in ( |3.8| ) and (3.9). 



(3.15) 



(3.16) 
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B. Average density of real roots 



On the real axi s the situation is rather different. There the imaginary parts of all fk and fk vanish due to the 
reality condition ( 1.1S| ). The rank of the matrix A is equal to 2 on the real axis so that the above calculations are 
invalid there. 

For small values of y = lm(z) the real and imaginary parts of the homogeneous polynomial F are 

Fi{x, y) = F(x, 0) + 0(y 2 ), F 2 (x, y) = yF'(x, 0) + <D(y 3 ), 



where F(x,0) is real^] and we use the prime to denote the derivative with respect to x. Introducing this into ( 1.16| ) 
— the definition of V 



V{a,a;z) = E j$(ai - F x (x, 0; o) + 0(y 2 ))6(a 2 - yF'(x, 0; o) + 0(y 3 )) 
x 8{a 1 - F'(x, 0; o) + 0(y 2 ))6 (a 2 - yF"(x, 0; o) + 0(y 3 )) 



(3.17) 



and comparing the arguments of the second and third deltas we see that the second one may be written as 8{a 2 
ya\ + 0(y 3 )). Setting the argument a 2 = 0, this becomes 

\a 1 \- 1 6( y ) + \y\- 1 6{a 1 ~o( y 2 )). 



Substituting the second t erm in the above sum into ( 3.1 7| ) would yield the s mall \y\ asymptotic behaviour of the 
density of complex roots ( 3.12 ) which we will obtain a little further on in [II C . 

The first term, which is zero off the real axis, generates the singular component of V{a\, 0, &i, a 2 ; x, y) 



^ S m 9 (ai,0, &i,a 2 ;z) = lr P (ai, &i; x)6(a 2 )6(y), 



(3.18) 



where Vq is the density for the joint distribution function of the real values of F (a;, 0) and F'(x, 0). 

In the Gaussian case Vo is also a Gaussian, determined by the 2x2 correlation matrix of the (real) values of F and 
F': 



A(x) 



(f(z),f(*)) (f(a0,f'0r)) 
(?(x),f(x)) (f'(x),f'(x)) 



(3.19) 



We use here the notations intr odu ced in the previous section for the real n- vector f (x), which has components fki(x, 0), 
and the bold round brackets (2.6) for the contraction over the Latin indices. 
Thus, in the Gaussian case: 



Vo(ai,&i;x) 



1 



27r x /det A(x) 



exp 



-l(a 1> a 1 )A-\x)(J 
2 \ai 



(3.20) 



Substituting the singular component P s i ng , with Vq given by ( |3.20 ), into (3.3) and performing the integral leads to 
Kac's formula Q for the average density of real roots in the homogeneous case 



£>o(x,j/) = S(y) 



(f(x),f(x))(f'(x),f'(x))-(f(x),f'(x))^ 



7r(f(x),f(x)) 



(3.21) 



Let us now obtain the singular component of the average density of roots for monic polynomials. Substituting 
Vsing, (3.18-3.2C), shifted by the deterministic part into (3.3) and performing the integration yields 



V^eTAR u{ J [0(x,O),0(x,O)] 

V (x,y) = 5(y) . H[w(x)\ exp <^ 

7r(f(x),f(x)j I 2(f(x),f(x)) 



(3.22) 



5 We remind the reader that the real and imaginary parts of a complex quantity are denoted by the indices 1, respectively 2. 
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where 



<>•)■■ <'M>(-y J + J^weii(w), (3.23) 



w(x) 



det A(x) 



and the error function is 



erf(s) 



(i(x),i(x)) 



#c,0) 



(3.24) 



~dt. 



The function H(w) is monotonically increasing. H(0) = 1 and H(w) w yjirj2w as io — > +oo. 

For homogeneous polynomials with non-central Gaussian distributions the density of real roots was obtained by 
Edelman and Kostlan [6]. Taking a singular limit of the covariance matrix in section 5 of [6] wi ll als o yield ( 3.22| ). 

Thus, in the Gaussian case the average density of zeros has a regular component given by ( [3.16 ) in terms of the 
2x2 matrices (/,/), (/,/) and (/,/)■ For the real type, (1.18), polynomials considered in this paper the average 
density of roots has also a singular component, localized on the real line, given by the generalized Kac type formula 
(fj§). 



C. Asymptotic behavior of the average density of complex roots near the real axis 



Let us now study the asymptotic behavior of the density of complex roots near the real axis. For small y = lm(z), 
the smallest eigenvalue of the matrix 



(/>/) 



(f,f) l/(f,f) 
y(f,f) y 2 (f',f) 



goes to zero as 



(f',f) — (fjf) /(f , f) y 2 . The argument of the exponential in the density of complex roots (3.16) 

goes to a finite limit; while the other factor in the numerator is 0(y 2 ). Thus, the average density of complex roots 
goes to zero as C(|y|): 



j(«^(x,o) *'l.*,oj)M*){p%%) , , 

V(x,y) = ^ { [(Qf,Qf)] + [Q<f>,Q<f>] ) + 0(\y 3 \), 

27T\/det A{x) L " > 

where the 2x2 matrix A(x) was defined above in ( |3.19| ) and the linear operator Q is defined by 

„ [(f , f) (f ', f ") - (f , f ') (f , f")] h! + [(f, f ') (f , f ") - (f , f ') (f ', f ")] h 



(3.25) 



Qh{x) = h" 



det A 



(3.26) 



We omitted all the arguments (x) of the functions appearing in the right hand side of ( 3.26|) . 

Thus, the real axis attracts the roots in its vicinity to the singular component located on it, depleting the density 
of roots in its neighborhood. 



D. High— and low-disorder limits 



With our definition of the moni c ra ndom polynomial (1.6) we may always consider that the random part of the 
polynomial (the coefficients Ck in fll.6| )) has expectation value equal to zero. Let us rescale the deterministic part <p 
of the polynomial to Tip, introducing a real parameter T > 0. The random polynomial is now 

P(z)=T(j>(z)+F(z). (3.27) 
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The form (3.27) allows us to interpolate from the monic to the homogeneous case. In the Gaussian case this is 
equivalent to considering the coefficients of the random polynomial as realizations of a Gaussian process with zero 
average and dispersion 



<y = r- 



(3.28) 



instead of (3.1). 

For small values of T — > 0, (the large rand omne ss limit, S — > oo) the average density of roots 2?, ( 3.16| ), approaches 
the de nsity for homogeneous polynomials, ( 3.12| ). In the large T limit, equivalent to small randomness, S — > 
in Eq.(3.28), the average density of roots concentrates near the roots of the deterministic part 4>{z) and decays 
exponentially away from them. Indeed, inspection of (3.16) shows that at fixed z = x + vy, which is not a zero of 0, 
the average density of roots T>, goes to zero exponentially when r — * oo. 

Let now zq = xq + iyo be a complex zero of <f> having multiplicity k > 1 : 



(z) = c(z-z ) k + O((z~z ) k+1 ), 



(3.29) 



and assume that the matrix A(zq) is non-singular. Then, to leading order in T 
p = \z — zq\, the average density of roots is 



-co and for small values of 



V(x,y) » Ak z T 2 p 



2 „2fe-2 e -i 



9 )r 2 p 2fc 



(3.30) 



where 



.4 



2 7 r v / A^' 



(3.31) 



% = arg(z - z ), 



B(0 ) = -c 2 [Ai + A 2 + |Ai - A 2 | cos(2fc0 o + Xo)] , 



Xa = arctan 



(f,f) n -(/,/), 



(3.32) 



(3.33) 



and Ai, A 2 are the (positive) eigenvalues of the 2x2 matrix (/, /). 

Now, it is readily seen that for T — > oo, (3.30) goes to k5 2 (z — z ) in the sense of distributions. For large but finite 
values of T the average density of roots is a Gaussian centered at zq for k = 1. 

For k > 1 if the eigenvalues of (/, /) are equal, Ai = A 2 , then the surface S = D(x,y) has an annular maximum 
for p 2k = (1 — fc _1 ) / BY 2 . If Ai ^ A 2 , the annular maximum splits into 2k individual maxima located at 2k9o + Xo — 
(2 M+ 1 )tt, M = 0, . . . ,2k— 1 and p 2k = (1 — k^ 1 ) / B m i n T 2 , where B m i n = c 2 \ m i n = mine B{6) is the minimal value 

of (gap- 

IV. ASYMPTOTIC DENSITY OF ROOTS FOR MONIC ALGEBRAIC POLYNOMIALS 

In this section we investigate the average density of roots for algebraic monic polynomials, 



P n (z) = z n + 6J2 



CkZ 



fc-i 



(4.1) 



fc=i 



in the large n limit. Here the paramete r 6 = > is the dispersion of the original Gaussian distribution which 
was transformed to the normal one , Eq.( [3.l| ), as mentioned in the preceding section. 
The roots of the polynomial (4.1) coincide with the roots of 



rp n (z) = rz n + J2ckz k -\ 



(4.2) 



k=l 



11 



so that we may use the rescaled <f>, (3.27), as in the preceding section. The large T (low disorder) asymptotic obtained 
there near complex roots of the deterministic part 4> cannot be used straightforwardly since in our case 4> has a highly 
degenerate real zero at the origin. 

The calculation of the 2x2 matrices (/, /), (/, /) and (/, /) in this case is presented in some detail in Appendix 
[A|. There we show that the matrix elements may be expressed in terms of the functions 

„ , , . 1 - r 2 cos 26 - r 2n cos 2nd + r 2n+2 cos(2n - 2)6 . . 

C n (r 2 ,29) = — 2 — — !s '-, 4.3 

1 — 2r z cos 26 + r 4 

2 r 2 sin 26 - r 2n sin 2nd + r 2 "+ 2 sin(2n - 2)6 

b ^ r ^ 20 >- l-2r 2 cos20 + r 4 ' (4 ' 4) 

and their derivatives with respect to r. Here 

z = x + iy = re lS . 

The matrix (/, /) is thus given by 

(f ft - 1 { C n{r 2 ,0) + C n (r 2 ,2ip) S n (r 2 ,2<p) \ 

V >:t >-2\ S n {r 2 ^) C n (r 2 ,0)-C n (r 2 ,2<p))- 

Its eigenvalues are C n (r 2 ,0) ± yj C 2 (r 2 , 28) + S 2 (r 2 , 26) . The smallest is positive for all 6 such that sin# ^ 0, as 
shown in 

In the Appendix we expressed the matrix eleme n ts of (/, /) and (/, /) in terms of derivatives of the matrix (/, /) 
with respect to r. We may now rewrite Eqs.( AlC - |All|) in matrix form using the orthogonal matrix 

v f cos 6 sin6>\ (4 6) 

w V — sin 9 cos 9 y ' 



(/=/) = 



(4.7) 



d(f,f) 



U T (6). 



(4.8) 



Here U T is the transposed of U. 

For examining the behaviour of the density of complex roots (3.16) in the vicinity of the unit circle, let us introduce 
a new, logarithmically rescaled variable 



e« =r 2 ". 



(4.9) 



Inspection of ( |3. 16 ) and comparison with (Q, Q, and (p^|) shows us that we need to estimate the 

matrices (/, /) and its inverse, the matrix 



U T (0)Tj 



I d(fj) 

2 dr 



(fjy\ 



and the second derivative with respect to r of the function C n (r 2 ,0). The functions C„(r 2 ,2<p) and S n (r 2 , 2<p) , (4.3) 
and (4.4), have rapid oscillations with 6. 



Let us study the quotient 

C„(r 2 , 2ip) (l - r 2 ) (l - r 2 cos 29 - r 2n cos 2n6» + r 2n+2 cos(2n - 2)6») 



Q 



C„(r 2 ,0) 



(1 - r 2 ") (1 - 2r 2 cos26» + r 4 ) 



for |£| <C The absolute value of the sum of terms proportional to r 2n in the numerator, which are fast oscillating 
with 6, does not exceed r 2n \J\ — 2r 2 cos 29 + r 4 . Noting that 1 — 2r 2 cos 29 + r 4 = (1 — r 2 ) 2 + (2r sin#) 2 , we obtain 



101 < 



1-r 2 



■ s r coth — , 

I J 2n|sin6»| 2 
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where we used the obvious inequalities 2ab < a 2 + b 2 and b < vV 
estimate can be obtained for the quotient S n (r 2 , 2(p)/C n (r 2 , 0). 
Thus, for 



b 2 to estimate the 9 dependent terms. A similar 



2n\ sinf 



■coth| < 1, 



(4.10) 



the matrix (/, /) is asymptotically proportional to the unit one times C n {r 2 , 0). Substituting the |(| < n asymptotic 
behavior of C„(r 2 , 0), ( |A20| ) we obtain that (/, /) is equal to n times a universal function of £, which does not depend 
on 9 and n: 



(/:/) 



n e^ 
2 ~ 



1 



1 + 0{t) 



(4.11) 



Noting that ( 4.10| ) implies £ <C n, we see that we may use (4.11) also for estimating the derivatives of (f,f) 



Let us first look at the homogeneous case, P n {z) =_XU=o c fe z " > corresponding to V — and recover the resul t 



obtained by Shepp and Vanderbei JhJ. Substituting ( 4.11 ) into (4/7), (L8), we obtain after a little algebra on ( 3.12| ) 

rf_d_ n 1 



P(eJ 



(4.12) 



Thus, for r< 1, (4.10), the expected density of complex roots for homogeneous (r = 0) Gaussian polynomials is 
asymptotically equal to n 2 times a simple universal symmetric function of £. For large values of |£| <C n the density 
has inverse power behavior, T> ~ £ -2 . 



In Fig. |l] we plot the renormalized average density of complex roots, n~~ D, for homogeneous real Gaussian 
polynomials of degree 99 as a function of £ = 2nln \z\ and 9 — arg(z). We see that at n = 100 the co ncor dance with 
the asymptotic formula (4.12) is good, excepting in the vicinity of the real axis where the condition ( 1.1C| ) is invalid. 

Performing the same substitutions on (3.16) and noting that for <f> = Tz 11 the 2-vector <j> satisfies 



4> = -U(0)</>, 
r 

we obtain after a little algebra a simple asymptotic formula, valid for t«1, for the average density of complex roots 
T> of monic algebraic polynomials: 



D(eJ 



n d 



i 



3« - 1 



exp 



r 2 



n 1 



(4.13) 



In the monic case the expected density of roots is also asymptotically equal t o n 2 times a universal function of £ 



and r jn. For nonzero values of r / 'y/n the large |£| <§; n asymptotic behavior of ( 4.13 ) is exponential decay for large 



positive £ and remains inverse-power for large negative values of £. This asymmetry becomes rather pronounced in 
the case of large values of the parameter r 2 /n, which was investigated numerically in Jl8j . 

Let us estimate n~ 1 J\f out (R), the fraction of the expected number of roots outside a disk of radius R centered at 
z = 0. For large n and | In i?| 2n the total number of real roots is O(logn) and since the sectors of angle C(7r/n) 
near the real axis contain a number of roots comparable to the error of the asymptotic relation (4.13), we may also 
use it there: 



Kut{R) 



1 



d£e' 



d9V(e^,9). 



Here £r = — \nR/(2n) and we may exploit the exponential decay of T> for large £ and replace e^'" by 1 under the 
integral if |£r| <C n. 

The fraction of roots outside a disk of radius R is thus asymptotically equal to 



-Mout{R) 

n 



1 



1 



2n \n(R) R 2n - 1 



exp 



2r 2 \n{R) \ 
R~ 2n ) ' 



1 



(4.14) 



For R — > this goes to n. 

In Fig. we plot the (re)normalized densi ty of r oots for n = 10 and N = T 2 = 100. It has nine sharp peaks in each 



of the half-planes Im(z) ^ 0. Our analysis n III D predicts the splitting of a complex n times degenerate zero of the 
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deterministic part into 2n maxima. In the present case we have only 2n — 2 peaks because the zero of the deterministic 
part is real and the other two peaks of the distribution are on the singular component. An interesting feature of this 
splitting is the fact that there are twice as many peaks as there are roots. Thus, on a typical realization of the random 
process we expect to find the n roots of the polynomial located near half of the positions of the maxima. 

At n — 30 and N = 1024 there is still some oscillation in the expected number of roots in an angular sector as can 
be see n in Fig. p[ For larger values of n, the density of complex roots approaches the 9 independent asymptotic form 
(4.13). In Fig. |f this convergence is shown for 2 | L <6'<-|at./V = r 2 = lOn. The number of extrema increases with 
n while their values become closer. The first several maxima near the real axis have a larger amplitude than those in 
the domain of universality. Since these are at a distance ~ n -1 from the real axis, the linear with |Im(z)| fall to zero 
of the density of roots may become rather steep as can be seen in Fig. 0. 



APPENDIX A: MATRIX ELEMENTS - ALGEBRAIC CASE 



The density of complex roots is given by Eq.( ft.l6 ) in terms of the 2x2 matrices (/,/), (/, /) and (/,/) defined 
in Section ||. Let us calculate them explicitly in the case when the polynomials are algebraic, i.e. when the basis 
functions are given by: 



f k (z) = z k -\ k = l,...,n. 

Taking z = re l6 , the components of the 2n- vectors / and / are 

Ai = r*- 1 cos(fc - 1)0, 
f k2 = r*- 1 sin(fc - 1)0, 
/ki = (fc-l)r fc - 2 cos(fc-2)0, 
fk2 = (k- l)r k ~ 2 sin(fc - 2)0 

Let us start with the 2x2 real symmetric matrix (/, /). 

n— 1 -. Ti—l 

(/, /) n = Y,r 2k cos 2 k9 = i £ r 2k (l + cos 2A;0) 



fc=0 



fe=0 



= -[C n (r 2 ,0) + C n (r 2 ,2<p)]. 
Here, we defined the function C n (x,x) as the real part of the sum 

n-l 

Y / x k e lkx = C n (x, X )+iS n (x, X ), 

1 — x cos x — x n cos nx + x n+1 cos(n — l)x 



k=0 

Cn(x,x) 



Sn(x,x) = 



1 — 2xcosx + x 2 
x sin x — x n sin nx + x n+1 sin(n — l)x 



1 — 2x cos x + x 2 

In a similar way, the other elements of the symmetric matrix (/, /) are given by 

(/,/) 12 = iS„(r 2 ,2 V ), 

{f,.f) 22 = \[Cn{r 2 ,0)-C n {r 2 ,2 V )l 

where the function S n (x 7 x) is defined in ( |A7| ). The determinant of the matrix (/, /) is 

det (/, /) = \ [C 2 n (r 2 , 0) - C 2 n (r 2 , 20) - S 2 (r 2 , 29)} 



1 — r 
1-r 2 



2n \ 2 



1 - 2r 2 " cos 2n0 + r 



in 



1 - 2r 2 cos 20 + r 4 



(Al) 

(A2) 
(A3) 



(A4) 

(A5) 

(A6) 
(A7) 



(A8) 



(A9) 



It is readily seen that the determinant is greater than zero if sin ^ i.e. off the real axis. 
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The sums appearing in the definitions of (/, /) and (/, /) are similar. The factors k and k 2 may be dealt with 

using the relation ^Z- = kx k ~ 1 . Thus, the elements of the matrices (/, /) and f.f , / ) ca n be expr essed in terms of 
derivatives with respect to the first argument of the functions C n and S n defined by ( |A6| ) and ( A7) above. 



(/>/)n 

(/>/)« 
(/>/)« 
(/./) 2a 



n-l 



fc=0 



[cos0(/, /) n +sin0(/,/) 21 ], 



^/cr 2 ^ 1 cos k9 cos(fc - 1)0 = J^fc? 

fc=0 

2<9r 

^|:[ co ^(/>/) 12 + sme(/ ! /) ; 
id_ r 

2d^ 

i_a r 

29^ 



,2fc-l 



[ cos cos 2 &;0 + sin 6* sin k9 cos fc0] 



'22 



(A10) 



cos0(/, f) 2l -sm0(f,f) n 
co S 0(f,f) 22 - S m6(f,f) 12 



n-l 

fc=0 
n-l 

- E fc 

fe=0 

i a a r 



„2 r 2fe-2 



cos 2 cos 2 k9 + 2 sin cos cos fc0 sin k9 + sin sin k9 



(A/) 12 
(A/) 22 



4r <9r dr 

_1__9 _9 

4r dr dr 



cos 2 9(f,f) n + 2 sin cos 0(/, /) 12 + sin 2 6(f, f) 2 
(cos 2 - sin 2 6) (/, /) u + sin cos ((/,/) 22 -(/,/) n 
cos 2 (/, /) 22 - 2 sin cos (/, /) 12 + sin 2 (/, /) u 



(All) 



(A12) 



After expressing all the relevant matrix elements in terms of the functions C n and S n , let us now consider their 
asymptotic behavior when n 3> 1. The sum in Eq.( |A5| ) and its derivatives with respect to x are truncated power 
series in x with power bounded coefficients. Their radius of conv erge nce is 1. For x < 1 the leading term of its 
large n asymptotic expansion is given by the n — > oo limit of Eq.(A5) yielding a result n-independent result with 
exponentially small corrections: 



Sn(x,x) 



1 — x cos x 
1 — 2xcosx + x 2 

a; sin% 
1 — 2x cos x + x 2 



0(x n ), 
0{x n ). 



(A13) 
(A14) 



If x > 1 the series the sum Eq.(A5) diverges as n — > oo. The leading terms for C n and S n are now exponentially 
large in absolute value and rapidly oscillating functions of x- The corrections are again rather small: 



C n (x,x) » x n 
S n (x,x) ~ x n 



x cos(n — l)x — cos nx 
1 — 2a; cos x + x 2 

a;sin(n — l)x — sinnx 
1 — 2x cos x + x 2 



+ 0(x- n ) 
f 0{x~ n ) 



(A15) 
(A16) 



Thus, in the large n limit the asymptotic behavior of the functions C n and S n c hang es dram atica lly in a narrow 
annulus of width C(n _1 ) near the unit circle, crossing over from the behavior Eq.(A13) to Eq.(A15) ( respectively 
from Eq.QA14D to Eq.flAl6|) ). 

For investigating the behavior in the crossover region, let us go to a logarithmic scale for the variable x: 



x = exp 



(!)• 



(A17) 
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For |£| -c nx, the leading asymptotic terms are 



i sin(n- |)x c g 







MY 



Sn(e»,x) 



y cos(n- \) x £ - 
cot Y e 

2 sin ^ 



(A18) 
(A19) 



Inspection of Eqs.( A18- AlS| ) shows that the corrections are small outside a neighborhood of the point x = 1, x = 0. 
For x = and |£| <C n, the leading term isQ 



C„ e-,0 



(A20) 



The asymptotic expansions given by Eqs.( [Al3| - |A16 ), (A18 - A20| ) may be now used to obtain the large n asymptotic 
behavior of the matrices (/,/), (/, /) and (/,/) in the domains \z n \ <C 1, \z n \ >■ 1 and the thin annulus which 
separates them. 



We will not write here the bulkier expressions which are valid for |^| < n without 
restrictions on x- 
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FIG. 1. The renormalized density of complex roots, n 2 D(z), for the homogeneous Gaus sian polynomial 
ciz" + C2Z 98 + . . . + cioo, n — 100, as a function of f = 2nln|z| and 8 = arg(z). (a) exact, Eq. (3.12); (b) asymptotic 



formula, Eq. (4.12); (c) contour plots: full line - exact, dotted line - asymptotic. 



FIG. 2. Renormalized density of complex roots for the polynomial N?z 10 + c±z 9 + . . . + cio, n = 10, and N = F 2 = 100, as 
function of £ and 8. 

FIG. 3. Average number of roots in angular sectors for n — 30 and N — 1024. Full line - histogram of the a verag es for 1000 



polynomials; dashed line - angular density of roots, J °° drrD(r, 8), obtained by numerical integration of Eq. ( 3.16) 



FIG. 4. Renormalized density of roots for N = T = lOn: (a) n = 100, exact; (b) n = 200, exact; (c) asymptotic formula 



(4.13) for N/n = 10. 



FIG. 5. Detail of the asymptotic behavior of the expected density of roots near the real axis for n — 30 and N = 1000. 



17 



